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Abstract: The transition from galactic to extragalactic cosmic rays is discussed. One of 
critical indications for transition is given by the Standard Model of Galactic cosmic rays, 
according to which the maximum energy of acceleration for iron nuclei is of order of _EpJ, ax ~ 1 x 
10 17 eV. At E > _EJ^ ax the spectrum is predicted to be very steep and thus the Standard Model 
favours the transition at energy not much higher than _BpJ, ax . As observations are concerned 
there are two signatures of transition: change of energy spectra and elongation rate (depth of 
shower maximum in the atmosphere X max as function of energy). Three models of transition 
are discussed: dip-based model, mixed composition model and ankle model. In the latter model 
the transition occurs at the observed spectral feature, ankle, which starts at E a ~ 1 X 10 19 eV 
and is characterised by change of mass compostion from galactic iron to extragalactic protons. 
In the dip model the transition occures at the second knee observed at energy (4 — 8) x 
10 17 eV and is characterised by change of mass composition from galactic iron to extragalactic 
protons. The mixed composition model describes transition at E ~ 3 x 10 18 eV with mass 
composition changing from galactic iron to extragactic mixed composition of different nuclei. 
These models are confronted with observational data on spectra and elongation rates from 
different experiments, including Auger. 

Introduction 

The Ultra High Energy Cosmic Ray (UHECR) 
has two most important problems. One of 
them is a presence of spectrum features pro- 
duced by propagation of UHECR particles 
through Cosmic Microwave Radiation (CMB) 
and the second is transition from galactic to 
extragalactic Cosmic Rays (CR). 

In the case of extragalactic protons two spec- 
tral signatures caused by interaction with 
CMB are predicted: Greisen-Zatsepin-Kuzmin 
(GZK) cutoff pQ and pair-production dip [2] . 

GZK cutoff is most spectacular prediction for 
UHECR, which status is still uncertain in 
present observations, though there are the in- 
dications to its presence. 

The pair-production dip is the spectral fea- 
ture originated due to electron-positron pair 
production by extragalactic protons interact- 
ing with CMB: p + jcmb — > P + e + + e~. 
Recently this feature has been studied in the 



works [31 El H] • The dip has been observed with 
very good statistical significance % 2 /d.o.f.^ 1 
by the Fly's Eye, Yakutsk, Akeno-AGASA and 
HiRes detectors, and with much worse statis- 
tical significance by Auger detector. 
The pair-production dip and GZK cutoff are 
signatures of protons. The confirmation of 
the shape of these features is the evidence 
for proton-dominated composition of primary 
CRs. For nuclei as primaries the shape of the 
dip and GZK cutoff are strongly modified. 

The different explanation of the dip has been 
proposed by Hill and Schramm [6]. They in- 
terpreted the dip observed in 1980s in terms of 
two-component model. The low energy compo- 
nent can be either galactic or produced by Lo- 
cal Supercluster. The similar model has been 
considered in [7|. The Hill-Schramm dip is 
widely used now for the explanation of the ob- 
served dip. 

From 1970s in the UHECR spectrum there was 
observed a flattening, which is called ankle. 
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Discovery of this feature at Haverah Park de- 
tector was interpreted as transition from the 
steep galactic component to more fiat extra- 
galactic one. The transition at ankle has been 
recently considered in [5j. 

In the dip model the transition is completed at 
the beginning of the dip at E « 1 x 10 18 eV. 
The ankle in this model appears as intrinsic 
part of the dip. Like in ankle model, the tran- 
sition occurs here also as intersection of flat 
extragalactic component (this flatness is espe- 
cially prominent in case of diffusive propaga- 
tion) with steep galactic spectrum. 

In the dip and ankle models the extragalactic 
component is assumed to be proton dominated, 
while the galactic component is most probably 
composed by iron nuclei. In the intermediate 
model, where transition occurs in the middle of 
the dip, the extragalactic CRs are assumed to 
have mixed composition [5]. 

In this paper all three above-mentioned models 
of transition are discussed. The logic of our 
discussion is as follows: we approach first the 
transition from the high energy end of galactic 
CRs, then we discuss the properties of UHECR 
relevant for transition problem and finally we 
describe the transition from properties of these 
two components. 

The end of galactic CRs 

With some disturbing small contradictions 
one may claim that at present we have the 
Standard Model for Galactic Cosmic Rays. 
It is based on Supernova Remnant (SNR) 
paradigm and includes four basic elements: (i) 
Supernova Remnants as the sources, (ii) SNR 
shock acceleration, (iii) Rigidity-dependent 
injection as mechanism providing the observed 
CR mass composition and (iv) Diffusive 
propagation of CRs in the galactic magnetic 
fields. 

(i) SNRs are able to provide the observed 
CR energy production in Galaxy, which 
can be found as Q ~ u> cr cM g /x CI [17], 
where us « 0.5 eV/cm 3 is the observed CR 
energy density, c is velocity of CR parti- 



cle, M g « 5 x 10 g is the total mass of 
galactic gas, and x cr « 7 g/cm 2 is the gram- 
mage traversed by CR before escaping from 
Galaxy. Using these numbers one obtains 
Q w 2 x 10 40 erg/s, which is less than 10% of 
energy release in the form of kinetic energy 
SNR ejecta per unit time. 

(ii) The great progress has been reached dur- 
ing last decade in the theory of acceleration. 
The cosmic ray streaming instability strongly 
amplifies the magnetic field upstream creat- 
ing highly turbulent field with strength up to 
SB ~ B ~ 10 -4 G [TU] (for recent works see 
[IT]). At each moment of the shock propaga- 
tion only particles accelerated to maximum en- 
ergy E max can escape outside. -E^iax reaches 
the highest value at the beginning of the Se- 
dov phase and then diminishes due to shock 
deceleration. The spectrum of escaping parti- 
cles has a narrow peak at energy E max (t) at 
each moment t, but the spectrum integrated 
over time has a classical E~ 2 shape with flat- 
tening at highest energies. This interesting re- 
sult has been recently obtained by Ptuskin and 
Zirakashvili [12"] . 

The maximum acceleration energy estimated 
in the Bohm regime of diffusion in the acceler- 
ation process is given by 

B / W \ 2 ^ 5 
£ max = 4 x 10 15 Z— — — ^ eV, 
10 4 G \n g /cm A J 

(1) 

where B is amplified magnetic field, W51 is the 
kinetic energy of the shell in units 10 51 erg, n g 
is upstream density of the gas and Z is charge 
number of accelerated nuclei. Thus for the pro- 
tons and iron nuclei the maximum energies are 

S max = 4 x io 15 J B_ 4 eV, 

££ ax = 1 x 10 17 B_ 4 eV. (2) 

E™ ax describes well the position of the proton 
knee and i?{JJ, ax predicts the position of iron 
knee. 

(iii) As the observations show, nuclei are sys- 
tematically more abundant in cosmic rays in 
comparison with interstellar medium in the so- 
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lar neighborhood [T3l HH [15] . The injection of 
particles in the regime of acceleration is respon- 
sible for it [T3J Q31 US] . It can be illustrated by 
simple consideration [TB] . 

A particle i from downstream (i — A, p) can 
cross the shock and thus to be injected in the 
regime of acceleration, if its Larmor radius 
tl{jp) > d, where d is the thickness of the shock 
front. Thus we readily obtain the relation be- 
tween nuclei and proton injection momenta 

P f ni = ZeBd/c=Zpf n .. (3) 

Eq. ([3]) results in < v? n ^, which provides 
the higher injection rate of nuclei. 

This conclusion can be reached also in more 
formal way. Consider flux of accelerated par- 
ticles i Ji{p) = Ki(p/p\ n ^)~ lg . Normalizing 
Ji(p) by condition 

/ Ji(p)dp = TJiHi, (4) 

c J pi . 

where rij is the density of gas i and r\i is a 
fraction of this density injected into accelera- 
tion process, we obtain for the ratio of fluxes 
of nuclei and protons in CRs 

Ja(p) _ zlg -i ?UnA _ ( 5 ) 
J p (p) ri p n p ' 

Thus, fraction of nuclei is enhanced by factor 
Ztja/Vp- For numerical calculations of CR 
nuclei abundances see [131 EH [Hj ■ 

(iv) CRs propagate in Galaxy diffusively, scat- 
tering off small-scale magnetic turbulence de- 
scribed as superposition of MHD waves with 
different amplitudes and random phases. This 
process is considered (see e.g. [17]) in the reso- 
nance approximation, when the giro-frequency 
of a particle is equal to a wave frequency in 
the system at rest with a motion of a particle 
along the average magnetic field. The mag- 
netic field is separated into average (constant) 
field B and fluctuating component B. In [T7] 
the parallel diffusion coefficient Dm (E) is cal- 
culated assuming D±(E) being much smaller 
than D\\(E) (see however the numerical sim- 
ulations [18] which does not support this as- 
sumption for the highest energies). 



The diffusion coefficient and its energy depen- 
dence is primarily determined by spectrum of 
turbulence w(k) which in most important cases 
is given in the power-law form w(E) oc k~ m , 
where k is a wave number. Then one has 

w(fc)ocfc~ m , D(E)<xE~ n , n = 2-m. (6) 

Thus, we obtain for the Kraichnan turbulence 
spectrum, which Landau and Lifshitz |19j con- 
sider theoretically preferable for MHD waves, 
to = 3/2 and D(E) oc E 1 / 2 - for the Kol- 
mogorov spectrum m = 5/3 and D(E) oc E 1 ^ 3 
and for diffusion in shock-dominated turbu- 
lence m = 2.0 and D(E) = const. In the cases 
when the turbulent magnetic component 6B 
is much larger than regular component Bq the 
Bohm diffusion is valid D(E) oc E; this case is 
in particular valid for acceleration on the shock 
fronts. 

Diffusive propagation is the only phenomenon, 
which imposes currently the problems for the 
Standard Model of Galactic CRs. The essence 
of this problem can be easily seen. 

Using the generation spectrum in the Galaxy, 
as that in acceleration Q gcn (E) oc E~ 2 , one 
obtains the diffuse spectrum: 

J(E) oc Q gcn (E)/D(E) oc £-( 2 +") . (7) 

Then from the observed spectrum J(E) oc 
E~ 2 - 7 one obtains D(E) oc E - 7 , which 
in principle results in too high anisotropy 
8(E) oc D(E) and too low traversed grammage 
x CI (E) oc 1/D(E) at high energy E. There are 
suggestions how these problems may be solved: 
the problem of small x CT (E) - by spallation in- 
side CR sources and reacceleration, the prob- 
lem of anisotropy - by local character of this 
phenomenon, and flat spectrum of helium - by 
acceleration in SNI remnant enriched by he- 
lium (see [2Ql[21] for discussion and references). 

The proton and nuclei spectra calculated re- 
cently by Berezhko and Volk [21] within Stan- 
dard Model are shown in Fig. [1] in comparison 
with observational data. The agreement with 
observations is quite good at low energies, and 
the knee is confirmed at i?kn ~ 3 x 10 15 eV in 
proton and all-particle spectra. The iron knee 
located at Ep e «8x 10 16 eV is most important 
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Figure 1: Fluxes and spectra calculated within Standard Model in [5T] for all particles, protons 
and nuclei are shown as function of kinetic energy efc. They are compared with data of CAPRICE, 
ATIC-2, JACEE and KASCADE. The position of the knees for all nuclei are given by e kn ~ 
3Z x 10 15 eV. The end of Galactic spectrum is given by iron knee ^ « 8 x 10 16 eV. At higher 
energies galactic spectrum becomes very steep. 



prediction of Standard Model. The spectra be- 
yond the knees are predicted to be very steep. 

Pair-production dip and GZK cut- 
off 

Being a quite faint feature, the e + e~- 
production dip is not seen well in the natu- 
rally presented spectrum log J(E) vs. logi?. 
The dip is more pronounced when analyzed in 
terms of the modification factor [H [3]. It is 
defined as a ratio of the spectrum J P (E) calcu- 
lated with all energy losses taken into account , 
and unmodified spectrum Jp nra (E), where only 
adiabatic energy losses are included. 

V (E) = J p (E)/J^(E) (8) 

The modification factor is presented in Fig. [2l 
If one includes only adiabatic energy losses, 



r]{E) — 1 according to definition (dash-dot 
line). If e + e~-production energy losses are 
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Figure 2: Modification factor 



additionally included one obtains dip, shown 
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Figure 3: The predicted pair-production dip in comparison with Akeno-AGASA, HiRes, Yakutsk 
and Auger data [22]. The first three experiments confirm dip with good x 2 /A.o.i. « 1.0 — 1.2, 
while the Auger data are characterized by larger x 2 /d.o.f. (see the text). 



in Fig. [2] by curve "e + e~". If to include 
the pion production, the GZK feature appears 
(curve "total"). The observable part of the 
dip extends from beginning of GZK cutoff at 
£«4x 10 19 eV down to E w 1 x 10 18 eV, 
where T] ~ 1. It has two flattenings: one at 
energy E a w 1 x 10 19 eV and the other at 
Eb rs 1 X 10 18 eV. The former automatically 
produces ankle (see Fig. [3]) and the latter pro- 
vides the intersection of flat extragalactic spec- 
trum at E < 1 x 10 18 eV with more steep galac- 
tic spectrum. 

The modification factor is less model depen- 
dent physical quantity than the spectrum. In 
particular it depends weakly on spectral index 
of generation spectrum 7 9 : In Fig.[2]the curves 
are plotted for 2.1 < 7 g < 3.0 with intervals 
A'jg = 0.1. The remarkable property of vis- 
ible dip in terms of modification factor is its 
universality. Modification factor rj(E) is given 



as dimensionless numbers for different energies 
and the curve remains the same when various 
physical phenomena are included in calcula- 
tions [4] : discreteness in the source distribution 
(the distance between sources may change from 
1 Mpc to 60 Mpc), different modes of propaga- 
tion (from rectilinear to diffusive), local over- 
density or deficit of the sources, large-scale in- 
homogeneities in distribution of sources, some 
regimes of cosmological source evolution (most 
notably those observed for AGN) and inter- 
action fluctuations. The only phenomenon 
which modifies dip noticeably is presence of 
more than 15% of nuclei in primary radiation. 
Therefore the proton dip in terms of modifi- 
cation factor is the universal spectral feature, 
determined mostly by interaction with CMB. 

The observed modification factor is given ac- 
cording to definition by the ratio of observed 
Jobs(E) to unmodified (J unm (E) cx _E~ 7 ") 
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Figure 4: The fluxes from Akeno-AGASA, HiRes 
tion. 

spectrum: ?y b s oc J h s (E)/ E~ ls , where 7 g 
is the exponent of the generation spectrum 
Qgen(Eg) oc Eg 7s in terms of initial proton en- 
ergies Eg. As Fig. [3] shows the pair production 
dip and beginning of GZK cutoff up to energy 
1 x 10 20 eV is reliably confirmed by all experi- 
mental data including AG AS A. As to AG AS A 
excess at E > 1 X 10 20 eV it can be explained 
by some other reasons, e.g. at some conditions 
by statistical fluctuations seen in MC of the 
work 23J. 

The comparison of the predicted dip with ob- 
servational data includes only two free pa- 
rameters: exponent of the power-law genera- 
tion spectrum j g (the best fit corresponds to 
7 9 = 2.6 — 2.7) and normalization constant to 
fit the e + e _ -production dip to the measured 
flux. The number of energy bins in the differ- 
ent experiments is 20 - 22. The fit is charac- 
terized by x 2 /d.o.f. = 1.0 - 1.2 for AGASA, 
HiRes and Yakutsk data. For the Auger data 
X 2 is good for hybrid data and very bad for 
surface detector data, mainly due to data in 
two lowest energy bins at 4.3 and 5.5 EeV. In 
Fig. [3] the hybrid spectrum shown by circles, 
and combined spectrum (surface detector data 
combined with fluorescent data at low ener- 
gies) shown by triangles are displayed. If to 
introduce the random energy errors SE/E in- 
side a bin (see section 'Discussion and Conclu- 
sions'), which is reasonable for the low energy 
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end of the surface detector measurements, \ 2 
is tremendously improved. The analysis will 
be presented somewhere else. 

One can see that at E < Ef, = 1 x 10 18 eV the 
experimental modification factor, as measured 
by Akeno and HiRes, exceeds the theoretical 
modification factor. Since by definition modifi- 
cation factor must be less than one, this excess 
signals the appearance of a new component of 
cosmic rays at E < Eb — 1 x 10 18 eV, and this 
component can be nothing but galactic cosmic 
rays. Thus, the transition from extragalactic 
to galactic cosmic rays, starts at energy E\>. 

The position and shape of the dip is robustly 
fixed by interaction with CMB and can be used 
for energy calibration of the detectors. 

The systematic errors in energy measurements 
are high, from 15% in AGASA to 22% in 
Auger. To calibrate each detector we shift the 
energies by factor A to reach minimum \ 2 m 
comparison with theoretical dip. We obtain 
these factors as A^ = 0.9, \y a — 0.75 and 
X m = 1-2 for AGASA, Yakutsk and HiRes 
detectors, respectively. After energy calibra- 
tion the fluxes given by AGASA, HiRes and 
Yakutsk detectors agree with each other in a 
very precise way (see Fig. |4} . The Auger flux 
is noticeably below the flux shown in Fig. [U 

In Fig. [5] we show in the left panel the com- 
parison of Auger data with that of AGASA, 
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Figure 5: Comparison of Auger data with AGASA, 
of Auger data and the energy shifted Auger data (. 
and Yakutsk data (right panel). 
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HiRes and Yakutsk. In the right panel we com- 
pare the Auger flux with the calibrated data 
of AGASA, HiRes and Yakutsk. We use also 
the energy-shifted Auger data (curve A = 1.2) 
with maximum shift allowed by systematic en- 
ergy errors of Auger. One can see that dis- 
agreement in fluxes survives. M. Teshima in 
his rapporteur talk [24] noticed that shift with 
A w 1.5 brings the data of Auger in agreement 
with the calibrated fluxes of AGASA, HiRes 
and Yakutsk. 

Three models of the transition 

In this section we describe three models of 
the transition from galactic to extragalactic 
CRs: ankle, dip and mixed composition mod- 
els. One feature is common for all three mod- 
els: they describe transition as intersection of 
steep galactic CR spectrum with more flat ex- 
tragalactic spectrum. One criterion which all 
models should respect is agreement with the 
Standard Model of Galactic CRs. The obser- 
vational data which has a power to confirm or 
reject each model include energy spectrum and 
mass composition. 

Ankle model 

This is a traditional model, based on the in- 
terpretation of the ankle as spectrum feature 
of the transition (see [8] for the recent works) . 



In fact this is most natural model, where tran- 
sition occurs because extragalactic component 
is very flat. This component is assumed to 
have pure proton composition with flat genera- 
tion spectrum oc E~ 2 valid for non-relativistic 
shock acceleration. Energy losses modify spec- 
trum insignificantly at E < 4 x 10 19 eV. The 
beginning of the ankle E a ~ 1 x 10 19 eV corre- 
sponds to equal fluxes of galactic and extra- 
galactic CRs at this energy. The transition 
at the ankle is illustrated by right panel in 
Fig. [BJ The curve "extr.p" presents the cal- 
culated extragalactic flux of protons and the 
dash-dot line gives the galactic CR spectrum. 
It is obtained by subtracting the extragalac- 
tic flux from the total observed flux following 
the procedure first suggested in [25]. The ob- 
served dip in the spectrum is explained not 
by pair-production dip, but by Hill-Schramm 
mechanism [6 . One must assume that galac- 
tic flux is presented by iron nuclei, and even in 
this case the ankle model contradicts the Stan- 
dard Model of Galactic CRs, since the half of 
the observed flux at E ~ 1 X 10 19 eV has the 
galactic origin. This model needs another com- 
ponent of galactic CRs with acceleration to en- 
ergy 100 times greater than maximum energy 
in the Standard Model. 

Another problem of this model is given by mea- 
sured elongation rate A max (_E), where A max is 
the depth of the atmosphere (in g/cm 2 ) where 
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Figure 6: Transition in the dip (left panel) and ankle (right panel) models. In both cases a solid line 
gives the calculated spectrum of extragalactic protons and a dashed line - spectrum of galactic iron. 
E tr is the energy of intersection of galactic and extragalactic spectra and Ep e gives the position 
of iron knee. Eb = 1 X 10 18 eV in the left panel is the energy where transition from galactic to 
extragalactic CRs is completed. 



a shower has maximum. In the right panel of 
Fig. E]X max (-E') calculated for the ankle model 
is plotted in comparison with elongation rates 
measured by different detectors. One can see 
that in energy range (1.5 — 5) x 10 18 eV there is 
great discrepancy between elongation rate cal- 
culated in all models with measurements of all 
detectors [25] . 
Dip model 

It is based on spectral confirmation of pair- 
production dip in energy range 1 x 10 18 — 4 x 
10 19 eV and beginning of GZK cutoff in energy 
range 4 x 10 19 — 1 x 10 20 eV. Since both of these 
features are signatures of protons, their ob- 
servational confirmation means the indication 
that mass composition is dominated by pro- 
tons. The shape of the dip allows admixture of 
nuclei not more than 10 - 15 %. The transition 
from galactic to extragalactic CRs is completed 
at E}, w 1 x 10 18 eV. The appearance of galac- 
tic CRs at E < Ef, can be seen from behavior 
of modification factor in AGASA and HiRes 
experiments below Eb (see Fig. [3]) and from 



flattening of calculated spectra for both recti- 
linear and diffusive propagation. The diffusive 
propagation makes flattening of the spectrum 
below Eb more pronounced [27[ [28] [29] . One 
can see this spectrum behavior for the case 
of the Bohm diffusion in Fig. [6] (left panel, 
curve "extr. p"); the apparent falling-down 
shape of this curve is caused by multiplication 
of the spectrum by E 2 - 5 . The intersection of 
this curve with the galactic spectrum, shown 
by dashed line, provides the transition from 
galactic to extragalactic CRs. The transition 
occurs at energy E tI ~ 5 x 10 17 eV. The galac- 
tic component is found by subtracting the cal- 
culated extragalactic proton flux from the ob- 
served flux, given by the KASCADE and HiRes 
data. Since the energy of transition E% T is 
close enough to the position of iron knee given 
by Eq. the dip model fits perfectly the 
Standard Model of Galactic CRs. The galac- 
tic spectrum below the iron knee is presented 
by iron nuclei, and thus transition takes place 
sharply between iron nuclei and protons. The 
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Figure 7: Elongation rate for the dip model (left panel) and ankle model (right panel). The 
calculated elongation rates are shown by the solid lines for QGSJET 30J model of interaction, 
by dashed lines for QGSJET-II and by dotted lines for SIBYLL [32 . The data points are 
measurements of Fly's Eye (stars), HiRes-Mia (squares), HiRes (circles) and Auger (triangles). 



feature in the observed spectrum, which corre- 
sponds to the transition in the dip model is 
the second knee at energy (4 — 8) x 10 17 eV 
as observed in different experiments. The as- 
sumed generation spectrum in the dip model 
has 7 9 w 2.6 — 2.7. Being extrapolated to 
E m i n i~ 1 GeV, such spectrum results in too 
high energy output of the sources. This prob- 
lem is naturally solved with an assumption 
that the actual source spectrum has a standard 
shape with j g = 2.0 for non-relativistic shocks 
or 7 g = 2.2 — 2.3 for relativistic ones. How- 
ever, the natural distribution of the sources 
over -E max |33j or luminosities [16] results in 
steepening of energy spectrum of generation 
rate Q(E g ) per unit volume of the universe to 
larger j g , starting from some energy. 

The prediction for elongation rate X max (E) is 
shown in Fig. (left panel). The characteris- 
tic feature of the dip model - sharp transition 
from galactic iron to extragalactic protons - 
results in steep increase of X max (E) with E 
below 1 x 10 9 GeV in contrast to the ankle 
model, where the increase of X max (E) is less 
steep, because of very flat proton spectrum at 
E < 1 X 10 18 eV. The observational data do 
not contradict the predicted steep increase of 
elongation rate below 1 x 10 18 eV. 

The left panel of Fig. [7] shows a reasonable 
agreement of the dip model with the bulk of 
experimental points in this figure, especially if 



one takes into account 20 — 25 g/cm 2 of sys- 
tematic error in all experiments. However, the 
detailed comparison of the dip prediction with 
the data of each experiment shows the different 
picture. While elongation rate predicted in the 
dip model agrees well with HiRes and HiRcs- 
Mia data, it does not agree with the Auger data 
especially with two highest energy points. 

The mixed composition model 
The main concept of the mixed composition 
model (see Allard et al from [8], [S], [34] . 
[35] ) is based on the argument that any ac- 
celeration mechanism operating in the gas in- 
volves the different nuclei in acceleration pro- 
cess and thus the primary flux must have mixed 
composition. For injection into process of 
acceleration the authors assume A-dependent 
regime, instead of rigidity-dependent one (J3J) 
in the Standard Model, and obtain J(E) cx 
^47 s -i£'-7 9 instead of Eq. Q valid for the 
Standard Model. It results in higher abun- 
dance of CRs by heavy elements in comparison 
with the Standard Model. In fact, as discussed 
in [16j . there are the reasonable regimes of in- 
jections when abundance of heavy nuclei is sup- 
pressed. The UHE extragalactic nuclei propa- 
gating through infra-red (IR) and CMB radia- 
tion are efficiently photo-disintegrated starting 
with energy E ~ 1 x 10 19 eV, while protons 
survive and therefore GZK feature is present 
in the mixed-composition model. At energy 
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Figure 8: The spectrum (left panel) and elongation rate (right panel) for the mixed model [53] 
with 7 g = 2.3, E a — 3 x 10 18 eV, cosmological source evolution (1 + z) 3 at z < 1.3 and a set 
of parameters Xi (see the text). At E > 3 x 10 19 eV the spectrum is strongly proton-dominated 
and is characterized by GZK cutoff (left panel). The mass composition evolves from almost pure 
iron composition at E w 3 x 10 17 eV to the lighter composition due to enrichment by protons 
and light nuclei of extragalactic origin. At energy E a = 3 x 10 18 eV the transition to pure 
extragalactic component is completed and chemical composition evolution proceeds further due 
to photo-disintegration of the nuclei. At energy E « 1.3 x 10 19 eV, seen in the plot, all nuclei 
are disappearing faster than before and composition becomes strongly proton-dominated at E > 
3 x 10 19 eV. 



below 1 x 10 17 eV the authors consider the 
mixed-composition spectrum |34j which is pro- 
portional to injection spectrum in the Galaxy: 

Qi(E) =x i Aj 3 - 1 KE-^, (9) 

where K is a normalization constant, i is a 

type of nuclei, free parameters, which 

describe the source chemical composition, and 
7 S is a spectral index, chosen to fit the data, 
with preferable values between 2.1 - 2.3, moti- 
vated by acceleration at the relativistic shocks. 
The cosmological evolution of the sources are 
included in calculations using factor (1 + z) m 
up to z m ax with different m, including m = 0, 
and different z max . 

With © taken as generation spectrum, the 
authors calculate the diffuse spectrum at 
higher energies propagating protons and nu- 
clei through IR and CMB radiation from the 
sources distributed uniformly in the universe. 
Using the calculated spectrum they fit the ob- 
served spectrum at energy higher than E a = 
3 x 10 18 eV, which is thus the energy where the 



pure extragalactic spectrum starts, i.e. transi- 
tion is completed. 

The galactic component is found by subtrac- 
tion of calculated extragalactic spectrum from 
the total observed spectrum. This procedure, 
adopted from [25], gives the spectrum be- 
low E a as observed and provides the smooth 
transition to calculated extragalactic spectrum 
at E a . Therefore, the part of the observed 
dip below E a is reproduced in this procedure 
phcnomenologically in contrast to the pair- 
production dip, which is accurately calculated. 

The calculated spectrum and mass composi- 
tion depends on parameters x\ in Eq. ([9]), 7 g , 
parameters of cosmological evolution m and 
Zmax, and therefore it is most flexible model 
among the three models at the discussion, 
which is able in particular to reproduce an 
arbitrary mass composition with some excep- 
tion. The robust prediction for spectrum and 
mass composition is related to energy range 
E > 1 x 10 19 eV, where the fraction of protons 
becomes large and steadily increasing, result- 
ing thus in the GZK feature and almost pure 
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proton composition, in contradiction with re- 
cent results of Auger. 

The spectra and mass composition predicted 
in one of the versions of the mixed model [34) 
are displayed in Fig. [8l The mass composition 
is in a good agreement with the selected data 
of Fly's Eye (only stereo), HiRes (only stereo) 
and HiRes-Mia, shown in the figure. The pre- 
dicted elongation rate has two break points, 
the first at E a = 3 x 10 18 eV, and the sec- 
ond at E « 1 x 10 19 eV. The first one occurs 
when transition to extragalactic CRs is com- 
pleted and evolution continues due to photo- 
disintegration of nuclei, first iron, then CNO 
and finally helium. At energy E > 1 x 10 19 eV 
(seen in the figure as E = 1.3 x 10 19 eV) 
all nuclei are destroying faster and at E > 
3 x 10 19 eV the composition becomes strongly 
proton dominated with GZK feature in the en- 
ergy spectrum. 

The first break point agrees with the Auger 
feature in elongation rate, but prediction of in- 
creasing A max at the second break point, i.e. 
at E > 1 x 10 19 eV, contradicts to decreasing 
of A max in the Auger data. 

Discussion and Conclusions 

The region of transition from galactic to ex- 
tragalactic CRs at energy between 1 x 10 17 — 
1 x 10 19 eV is the key energy range for under- 
standing the origin of CRs. At low energy part 
it includes the high energy end of galactic CRs. 
The information on maximum energy of accel- 
eration, chemical composition and propagation 
in Galaxy at these energies will clarify the to- 
tal picture of origin at lower energies. The low 
energy part of UHECRs is important for under- 
standing of origin of UHECRs and their prop- 
agation in extragalactic magnetic fields. The 
transition from galactic to extragalactic CRs 
is the central issue of this energy region. 

There are two detectors which cover partially 
the above-mentioned region: KASCADE- 
GRANDE [36] and TALE [37]. There are also 
the proposals to extend the observations of 
Auger to energy E ~ 1 x 10 17 eV (see e.g. [55]). 
The Auger detector has great potential to ex- 



plore this region, building more dense part of 
the detector covered with fluorescent, scintilla- 
tor and muon detectors. 

The basic information which can be obtained 
includes precise measurement of energy spectra 
and mass composition (there is little hope to 
detect anisotropy in this energy region, though 
in some models the galactic sources can be ob- 
served in protons with energy E < 10 18 eV 

At present we have the sufficiently good data 
on spectra and mass composition at energy 
range 1 x 10 18 - 4 x 10 19 eV. The spectra 
are measured with high statistics (especially in 
case of the Auger detector), but problem is the 
accuracy of energy determination. From quite 
disappointing Fig. [5] (left panel) one concludes 
that scales of energy determination is quite 
different in all detectors. Energy calibration 
with help of the pair-production dip suggests 
that energy measured by scintillator detectors 
is systematically higher than that by the flu- 
orescence detectors and it gives a reasonable 
recipe of increasing energies given by fluores- 
cent method and decreasing it for the scintilla- 
tion method. In this case the curves 'Yakutsk' 
and 'Akeno-AGASA' in Fig. [5] go down and 
'HiRes' and 'Auger' - up. For HiRes, AGASA 
and Yakutsk the method of calibration with 
help of dip works successfully (see Fig. 01 with 
energy shift within the allowed systematic er- 
rors, but for Auger it requires the shift by fac- 
tor 2 greater than systematic error. 

The pair-production proton dip in terms of 
modification factor is an excellent tool to mea- 
sure spectrum shape independently of abso- 
lute flux. From Fig. [3] one sees the excellent 
agreement of the theoretical dip with data of 
AGASA, HiRes and Yakutsk. By the stan- 
dards of cosmic-ray physics the agreement with 
Auger data is also good, but \ 2 f° r compari- 
son with SD data is very large. This is a re- 
sult of very big statistics in the surface detec- 
tors at lowest energies E > 4.5 x 10 18 eV. In 
the lowest energy bin at E — 4.5 x 10 18 eV 
there are 4128 events and the error in determi- 
nation of flux provided mostly by this statis- 
tics is SJ/J — 0.024. The theoretical value of 
modification factor at this energy is only 14% 
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higher than experimental value, but owing to 
very small S J/ J, the contribution of this bin to 
X 2 is 99.27 ! Most probably the other sources 
of errors should be included in the bins with 
small SJ/J, and a possible source of this error 
is the energy errors which are changing ran- 
domly inside a bin. These could be statistical 
errors and energy-dependent part of system- 
atic errors. Assuming that number of events 
are distributed in a bin as N(E) = KE~ 7 one 
obtains SJ/J = y(8E/E) r , where (SE/E) r is 
the random energy error inside the bin. The es- 
timated value 5 J/ J is much larger than what 
obtained in Auger analysis for all reasonable 
values of (SE/E) r and 7. More generally, ac- 
cording to Markus Roth's remark, \ 2 analysis 
is not adequate for the cases of small SJ/J and 
large (SE/E). At this stage of analysis we do 
not consider Fig. [3] as contradiction with Auger 
data. 

Coming to the transition from galactic to ex- 
tragalactic CRs, we emphasize that at present 
there are only two experimental methods to 
study it: measuring the spectrum and mass 
composition. The transition will be clearly 
seen if spectrum of iron nuclei and that of 
protons are measured separately (see Fig. [5]), 
but even without this ideal possibility the to- 
tal spectrum has signatures of transition in 
the form of the spectral features - second knee 
in case of the dip model and ankle in case 
of the ankle model. The spectrum can be 
measured nowadays with high accuracy and 
its shape contains the information about mass 
composition, which is the other characteris- 
tic of the transition. The pair-production 
dip with its specific shape is a signature of 
proton-dominated composition (nuclei contri- 
bution should be not more than 10 -15 % [5]) 
and its observational confirmation is an argu- 
ment not weaker than that due to A max mea- 
surement (we remind that only two free pa- 
rameters are involved in describing about 20 
energy bins in each experiment). 

The mass composition gives another way to 
test the transition. The best method at 
present is given by measuring of elongation 
rate X max (E). Unfortunately this method has 



many uncertainties, including those in value of 
fluorescent yield, absorption of UV light in the 
atmosphere and uncertainties in the models of 
interactions, needed to convert the tested mass 
composition into A" max . The systematic errors 
in measuring X max can be as large 30 g/cm 2 to 
be compared with difference about 100 g/cm 2 
between X max for protons and iron. The better 
sensitivity for distinguishing different nuclei is 
given by distribution over A max |26j . 

There are three models of the transition: ankle, 
dip and mixed-composition model. They differ 
most notably by the energy of transition (an- 
kle: E ~lx 10 19 eV, dip: B«lx 10 18 eV and 
mixed composition model E « 3 x 10 18 eV), 
and by mass composition of extragalactic com- 
ponent (protons - for the ankle model, proton- 
dominated - for the dip model and mixed com- 
position - for the third model). 

The ankle model contradicts the Standard 
Model of Galactic CRs (energy where galac- 
tic flux is half of that observed is two orders 
of magnitude higher than energy of iron knee) 
and severely disagrees with A max measured in 
all experiments at (1.5 — 5) x 10 18 eV. 

The dip model is based on well confirmed sig- 
nature of proton interaction with CMB - pair- 
production dip. The two other models must 
assume that agreement of pair-production dip 
with data is accidental and the observed dip 
is produced by two components, galactic and 
extragalactic. The dip model assumes the iron- 
dominated galactic flux below 5 x 10 17 eV 
and proton-dominated extragalactic flux above 
1 x 10 18 eV. This mass composition is con- 
firmed by HiRes and HiRes-Mia data for elon- 
gation rate. It does not contradict the bulk of 
all data on X max , but contradicts X max mea- 
sured by Auger, especially the highest energy 
points. The generation spectrum in this model 
is E~ 2 or E~ 2 2 as needed by shock accelera- 
tion with a steepening to 7 9 = 2.7 due to distri- 
bution of sources over maximum energy of ac- 
celeration of source luminosities. The proton- 
dominated composition can be produced in 
some models of injection to the shock accel- 
eration. 

The mixed composition model assumes mixed 
composition generation spectrum for extra- 
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galactic component with generation index 2.1 - 
2.3. It has many free parameters, most notably 
ones describing the mass composition of the 
generation spectrum, and thus it can in prin- 
ciple explain any observed mass composition. 
However, this model has a robust prediction 
at energy E > 3 x 10 19 eV: proton-dominated 
composition and the GZK feature. As far as 
Auger elongation rate is concerned, the mixed 
composition model explains well the break in 
elongation rate at 2 x 10 18 eV and contradicts 
the two Auger points at E > 2 x 10 19 eV. The 
energy where transition to extragalactic CRs is 
completed in most versions of this model equals 
E«3x 10 18 eV. Much better quality of data 
on X max is needed to distinguish the dip and 
mixed-composition models by A max measure- 
ments. Probably it is possible to do using A max 
distribution [2"6] , 

We will comment now on agreement of the 
transition models with the measured galactic 
spectrum. For all three models it is reached by 
the formal subtraction procedure: the galac- 
tic spectrum is found as difference between 
measured total spectrum and calculated ex- 
tragalactic spectrum. But the galactic spec- 
trum calculated in the Standard Model at E > 
1 x 10 17 eV is very steep and, as was demon- 
strated in [39 , for diffusive model of propaga- 
tion all three models contradict the calculated 
galactic spectrum, the dip model to the less 
extent. Strictly speaking this contradiction is 
produced by exponential cutoff in the acceler- 
ation spectrum at E > 

The most consistent conclusions on nature of 
observed UHECRs are obtained at present 
by HiRes detector: it has confirmed the 
pair-production dip and thus proton-dominant 
composition at 1 x 10 18 - 4 x 10 19 eV, 
the X max measurements agree with proton- 
dominant composition at E > 1 x 10 18 eV , and 
E x j2 measurement confirms that steepening of 
the spectrum observed at E > 4 x 10 19 eV is 
really the GZK cutoff. Therefore, according to 
these data CRs observed at E > 1 x 10 18 eV are 
extragalactic protons exhibiting two signatures 
of interaction with CMB: pair-production dip 
and GZK feature. 
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